DIFFUSION DUE TO THE BEAM-BEAM INTERACTION AND FLUCTUATING FIELDS IN 

HADRON COLLIDERS 



D 

in 



> 

O 

o 

O 

\o 

Oh 
I 

o 
o 



Tanaji Sen 

DESY, 85 Notkestrasse, 22607 Hamburg, Germany 

James A. Ellison 
University of New Mexico, Albuquerque, NM 87131 

Random fluctuations in the tune, beam offsets and beam size in the presence of the beam-beam interaction are shown 
to lead to signiflcant particle diffusion and emittance growth in hadron colliders. We find that far from resonances high 
frequency noise causes the most diffusion while near resonances low frequency noise is responsible for the large emittance 
growth observed. Comparison of different fluctuations shows that offset fluctuations between the beams causes the largest 

diffusion for particles in the beam core. 



Emittance growth due to the non-linear beam-beam in- 
teraction is a major concern at all hadron colliders. The 
tunes (i.e. the rotation numbers of the transverse beta- 
tron oscillations) are always chosen to avoid low order res- 
onances and the overlap of high order resonances does not 
create chaotic regions large enough to lead to significant 
amplitude growth. The dynamics are qualitatively differ- 
ent however when the time-dependence, both deterministic 
and stochastic, of parameters is included. The effects of 
.deterministic tune modulation are well studied and remov- 
ing modulation lines from the betatron spectrum reduces 
particle loss from the tails of the beam Q . Random fluc- 
tuations in the tune, closed orbit and beam size are also 
present in accelerators. Qualitative arguments and nu- 
merical simulations |^ have shown that tune fluctuations 
lead to emittance growth especially for tunes close to a 
■resonance . Our aim is to provide a quantitative theory to 
explain these observations. 

We consider collisions of a proton beam with an op- 
posing beam composed of either leptons as at HERA or 
hadrons as at the Fermilab Tevatron and the proposed LHC 
at CERN. For one dimensional motion, the beam-beam 
potential seen by a proton, assuming a Gaussian charge 
distribution of the opposing beam, is given by V{x) = 
Cl^ dq \\ ~ e^^'/(2<T^p+'?)]/(2a2p q). The constant is 
C = Nb^opTp/^p where N},^op is the number of particles per 
bunch in the opposing beam, rp is the classical radius of 
the proton, 7p is the relativistic kinematic factor for the 
protons, and Uop is the sigma of the size of the opposing 
beam. Transforming to action-angle coordinates (J, ?/') via 
X ~ \/2J/3* cos tA, x' — ~^2J/ (3* sin V' where fi* is the 
(3 function at the interaction point, we obtain the Fourier 
expansion of the potential 



U{x) = C[Fo (a) 2 ^ Fk (a) cos 2kil)] 
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where a = /3*J/(2(T^ ) is a dimensionless amplitude. The 



where the Ik are modified Bessel functions and 2F2 is a 
generalized hypergeometric function. Including the linear 
motion and the beam-beam interaction, the Hamiltonian is 
H = v^J + U{J,ip)5p{9) where Sp{9) is the periodic delta 
function with period 2tt/Nip, Njp being the number of 
interaction points and 0, the "time" variable, advances by 
27r per turn. Integrating the equations of motion over one 
turn leads to the one-turn beam-beam map: A-)/) = 27ri^'^ -I- 
dU/dJ and A J = -dU/d^^. 

First we consider the diffusion in amplitude due to ran- 
dom fluctuations in the tune. In the extreme case when the 
tune is completely random on the interval [0, f], the beam- 
beam kicks occur at uncorrelated phases and the emittance 
grows as in a random walk process. Usually the random 
contribution to the tune is quite small, of the order of 0.001 
at the most, but this is sufficient to affect the long time dy- 
namics. The sources of tune fluctuation include power sup- 
ply noise in quadrupoles, closed orbit fluctuations through 
the non-linear magnets and mechanical vibrations of the 
non-linear magnets. In addition, both intra-beam scatter- 
ing due to the Coulomb force and RF noise lead to fluctuat- 
ing particle momenta. This in turn leads to a tune fluctua- 
tion via the machine chromaticity. We model the tune fluc- 
tuation by an additional term At/j^ in the total phase. As- 
suming that the random contribution is small, we can write 
the change in action at turn m due to this fluctuating phase 
alone as AJr{m) = [dAJ{m)/dij] A^r{m) + 0{Ailjf). The 
unperturbed total phase at turn m and action J is ip{m) = 
2mT:v{J) -f tpQ^ where V'o is the initial phase. When the 
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tunc is far from a resonance, the linear action J is con- 
served after averaging. This allows us to assume that 
AJr(m) <^ t^(0) so that in the sum over turns we can re- 
place J(m) by J(0). We assume that the random process 
is stationary so that the random phase correlation func- 
tion is of the form (AV'r(O^V'r(n + 0) = ^Tr'^^i'r K,y{n), 
where the average is over many realizations of the noise 
process, Av^ is the amplitude of the tune fluctuations and 
K^{—n) — K^{n). The diffusion coefficient defined as 
D^iJ) = liiUN^ooiiJiN) - J(0)]2) /N is found, by extract- 
ing the dominant terms, to be 
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We observe that only the tune noise at even harmonics 
of the betatron tune leads to a diffusion in the action. 
Hence far from resonances, high frequency components of 
the tune fluctuation spectrum cause diffusion. A natu- 
ral choice to model the tune fluctuations is the Ornstein- 
Uhlenbeck (OU) process because it is the only Gaussian 
stationary Markov process and the spectral density S'(w) 
(related to by the cosine transform) decays as w"^ 
which is in reasonable agreement with measured noise den- 
sities. For the discrete time OU process with correlation 
time Tc, the correlation function is K^{n) = (1— 1/tc)"/[1— 
l/(2rc)]. The spectral density drops to roughly half its 
maximum value at a frequency /1/2 — frev/i'^'^Tc)- The 
revolution frequency frev is 47.3kHz for HERA. Substitut- 
ing this form for K^^ leads to 



L'^(J) = 128 



1 - 1/(2tc) ^ cosh9 - cos47^fc^ 



(5) 



where Q = — ln(l — I/tc). The main amplitude dependence 
of D^{J) is contained in the Fourier coefficients. The dom- 
inant coefficient Fx has the expansion Fi = '^^i{—^y^^ 
(a/2)-'{(2j)!/[(j-hl)!(j!)2]}. At small amplitudes therefore, 
D^{J) ~ J^. At large amplitudes the beam-beam force 
vanishes and all Fourier coefficients Fk go to constant val- 
ues . Hence the diffusion coefficient D^{J) increases mono- 
tonically in the core of the beam and levels off in the tails of 
the beam. Measurements of beam loss rates at HERA |^] 
are consistent with a constant diffusion in the tails of the 
beam. The mean escape time to an action Jf,, calculated as 
Tesc{J) = J dJ/D^{J), goes to infinity logarithmically 
at small amplitudes as Tesc{J) ~ ^^{Jb/ J)-, while at large 
amplitudes it goes to zero as Tesc{J) ~ {Jh — J^)- 

We have compared the above analysis with a numeri- 
cal calculation. An initial distribution of f 000 particles is 
placed at 100 different amplitudes with 10 particles at each 
amplitude distributed uniformly in phase. The particles 
are tracked for 10^ turns or more (the number increasing 
with the noise correlation time) using the beam-beam map 
with the random contribution to the tune at each turn 
determined by an OU process. The diffusion coefficients 
at each amplitude are averaged over the ten phases and 
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Figure 1: The diffusion coefhcient D^{J) calculated the- 
oretically from Eq. (5) (dashed line) compared with the 
values obtained from the simulation. Parameter values: 



=0.291, Ai^r 
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over ten noise realizations. The simulations were done for 
three correlation times: Tc = 1.1, 10, 100. Figure 1 shows a 
comparison of the diffusion coefficient obtained from Equa- 
tion (H) with that calculated numerically at Tc = 10. The 
level of agreement is nearly the same at other correlation 
times. Both the analysis and the numerical results show 
that D^{J) ~ for large Tc. High frequency noise (above 
1 kHz say) is strongly attenuated within the interior of the 
beam pipes by the impedance of the magnets and other 
devices in a storage ring so realistic values of the noise 
correlation time for HERA are Tc > 10. 

The diffusion in amplitude causes emittance growth 
over the period of stored beam - typically 24 to 30 hours 
for the proton beam at HERA. The time evolution of the 
density and the emittance can be followed by solving the 
Fokker-Planck equation. Assuming that the diffusion in 
action is a Markov process and the drift coefficient is half 
the derivative of the diffusion coefficient (as is usual for a 
Hamiltonian system) j^, the Fokker-Planck equation for 
the density is 



dp 
dt 
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2dJ 



D{J) 



dp 
dJ 



(6) 



We integrate this one dimensional Fokker-Planck equation 
by the method of lines An absorbing boundary is 

placed at an action Jt, corresponding to the position of the 
beam-pipe. The density at the origin does not change since 
the diffusion coefficient and its derivative vanish there. 
The evolution of the average action is then found from 
{Jit)) ^ Jp{J,t)dJ/ f^" p{J,t)dJ. Figure 2 shows the 
evolution of the average action over nearly 30 hours for 
three correlation times. For Tc = 1.1, (J) grows the most 
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Figure 2: Relative growth in the average action (J) of a 
proton beam over a storage time of 30 hours due to tune 
fluctuations at three noise correlation times. Parameter 
values: ;/°=0.291, Az/^ = IQ-*. 

rapidly as expected, then decreases as particles are lost at 
the boundary. Diffusion is slower for the other two correla- 
tion times and particles are not lost at the boundary so (J) 
grows almost linearly with time. These calculations show 
that even in this one dimensional model, tune fluctuation 
can cause the emittance to grow by 10 to 70% over the 
storage time in a proton machine. 

Emittance growth due to fluctuations is significantly 
enhanced near a resonance. As the tunes approach a res- 
onance, the resonance islands increase in width and the 
resonant amplitudes move further out. Finally at the reso- 
nance tunes the fixed points have moved off to infinity be- 
cause the beam-beam tune shift is largest at the origin and 
vanishes at infinity. The phase portraits are strongly tune 
dependent. For example, with v — 0.25, the phase space 
portraits are diamond shaped close to the origin and at 
large amplitudes they are four-armed stars with long arms 
along the four axes. A level curve in the phase portrait 
can be labelled by an energy to be defined below. Particle 
motion with external fluctuations has two aspects: motion 
on a level curve and diffusive motion between level curves. 
Near resonance a particle initially close to the origin may, 
by diffusion, be eventually transported to a resonance is- 
land where it experiences a large jump in amplitude. Ex- 
actly on resonance, e.g. at v — 0.25, a particle may, after 
a long time, diffuse on to a star shaped curve which subse- 
quently leads to a very large amplitude excursion. Particles 
are most likely to leave a given level curve in the vicinity of 
the turning points since most time is spent in this neigh- 
bouhood. The projection of diffusive motion along the gra- 
dient to the level curves then transports particles to larger 



amplitudes. 

To analyze the diffusive motion we observe that without 
noise and even after averaging the linear invariant J is not 
conserved near resonance. Instead, after averaging over 
the fast varying phases, a time-independent Hamiltonian is 
obtained which describes motion close to the 2A;th integer 
resonance (tune i>2k = Integer/2fc), 

oc 

H = SJ + A[Fo + 2 '^F„ik cos2mk<j>] . (7) 

m— 1 

where S ^ ly'^ ~ h'2k, the difference from the resonance tune, 
(j) = -0 — i/2fc6', the slowly varying phase, and A = C/{2tt). J 
oscillates periodically betweeen two limits Jmin and Jmax 
which are determined by the transverse energy E = H. 
Dropping all Fmk, m > 1, the period on a curve labelled by 
EisTE = (l/2^fc) dJ/y^AA^F'^ -{E-SJ- AF^)^ 
In contrast to the betatron tune v, this tunc ve = ^/Te 
is very small, typically of the order of 0.001 for fc = 2 
and S ~ 10~^. i'e increases with the resonance order 
2k. The tune fluctuations Avr cause the resonant am- 
plitude and the island widths to also fluctuate and the 
energy to diffuse. The rate of change of H is found to be 
dH/dO — '-{dJ/d6)uAi/r, the subscript u denotes the un- 
perturbed rate of change. After integrating over a turn, 
the total change in the Hamiltonian at turn N is H{N) — 
H{0) = -J2l=i[Jirn) - J{m - l)]Aurim). Expanding J 
in a Fourier series: J{m) = X^^o cos{2tt veJ in + and 
using the stationarity of {Ai'r{l)Ai>r{l + we obtain for 
the diffusion of the energy 

D (E) = (A^.)^ - i3|(l-cos2^j^^)sinh9 
' l-l/{2Tc)j-[ coshQ ~ cos2-KjvE 

The Fourier amplitudes Bj grow as the tune approaches 
the resonant value from below. The diffusion at a given 
tune increases smoothly moving out from the origin, jumps 
when the particle is on the largest of the resonant islands, 
decreases to zero at the stable fixed point, increases back 
to the value on the largest island and stays nearly constant 
thereafter. The noise frequencies which contribute to the 
diffusion in energy are the harmonics of the low frequency 
t^Efrev The topology of the phase space orbits and the fact 
that noise of comparatively low frequencies has the dom- 
inant contribution to the diffusion in energy explains the 
large growth in emittance due to noise in the neighbour- 
hood of a resonance. 

For the ID beam-beam interaction, diffusion is signif- 
icantly enhanced only close to low order resonances, but 
for the 2D interaction even relatively high order resonances 
e.g. 14th order can lead to large emittance growth, as was 
observed in J3] . A detailed study of the 2D case will appear 
separately 

Next we consider fluctuations of the offset between the 
beams at the IP. This has two effects. The position of 
the maximum of the beam-beam force fluctuates so more 
particles in the proton beam will be subjected to a larger 
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force. It also destroys the symmetry of the beam-beam 
force and can excite odd order resonances. This had been 
studied earher by Stupakov |^ assuming the potential for 
a flat beam. Our results below assume the general form 
of the potential given in Eq. (|l|). We assume that the 
offset fluctuation d,.(rn) at turn m is small and write it 
as dr{m) — A(irx(™)o'op where Adr is the dimensionless 
amplitude of the offset and is a random variable of 

zero mean and unit variance. Calculation of the diffusion 
coefficient far from resonances yields 

-. oc 

DoffiJ) = -(Ca„pAd.)2^(2fc + I)2G2(a) 

oo 

X ^^Koffin) cos2Tr{2k + l)vn (9) 

n— — oo 

The correlation function is Koff = {xiOxi^ + 0)- ^k, the 
Fourier coefficients of the beam-beam force, are given by 
Gk = MK+i + K]/^op + [{k + l)Fk+i - kFk]/{V^aop). 
Notice here that the odd harmonics of the betatron tune 
contribute to the diffusion in action. 

We now discuss the diffusion due to beam size fluctua- 
tion of the opposing beam. The location of the maximum 
of the beam-beam force is oc Uop while the maximum itself 
is (X a~p. Consequently when Uop fluctuates, protons in 
a larger range of amplitudes will be subject to the maxi- 
mum of the force - as with offset fluctuations. A study of 
this for the flat beam potential was reported recently in 
1^. We find that the diffusion coefficient for the general 
beam-beam potential is 

oo oo 

D„^^{J)^32{CAar a)2^[fci^^]2 ^ K„^^ cos 'iirkvu 

k—1 n— — oo 

(10) 

The size fluctuation at turn n is A<Tr'n{n)Uop where Acr^. is 
dimensionless, rj{n) is a random variable of mean zero and 
unit variance, and Ka^^{n) = {ri{l)rj{n + I)). 

Figure 3 compares the analytical diffusion coefficients 
from the three fluctuating phenomena considered here. The 
parameter values are shown in the caption. We find that 
the diffusion due to offset fluctuations is largest for ampli- 
tudes less than 2ap. At greater amplitudes, diffusion due to 
beam size and offset fluctuations, both of which directly af- 
fect the amplitude of the beam-beam kick, are of the same 
order of magnitude [lO|] . Diffusion due to tune fluctuations 
is the smallest at all amplitudes because it affects only the 
phase at which the particle is kicked. Nevertheless, the 
sources of tune fluctuations are difficult to eliminate and 
more numerous than for the other fluctuations. 

To summarise, we have reported three main results in 
this letter. Far from low order resonances, high frequency 
tune fluctuations cause larger growth of particle amplitudes 
than low frequency fluctuations. These high frequency fluc- 
tuations can cause the emittance to nearly double over the 
storage time of a day. Near resonances, low frequency fluc- 
tuations are resonant with the motion of the linear invari- 
ant and these lead to the largest diffusion in the energy 
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Figure 3: Comparison of the diffusion coefficients in J due 
to tune fluctuations, offset fluctuations and beam size fluc- 
tuations given by Equations (5), (9) and (10) respectively. 
Parameter values: i^°=0.291, = 10.0, Aur = 10"'*, 
Adr = 0.01 = Aur- 

which subsequently leads to signiflcant emittance growth. 
Comparing different fluctuations in the off-resonance case, 
we have found that for reasonable values of the fluctuat- 
ing amplitudes, offset fluctuations at the interaction points 
cause the largest diffusion at small amplitudes while at 
large amplitudes, fluctuations in the size of the opposing 
beam have a comparable effect as the offset fluctuations. 

We thank A. Bazzani, H. Mais and F. Willeke for very 
fruitful discussions and support. 
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